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Abstract 

The aim of this paper is to discuss the role of robustness and diversity in population dynamics in particular 
to some properties of the multi-step from healthy tissue to fully malignant tumours. Recent evidence 
shows that diversity within the cell population of a neoplasm, a pre-tumoural lession that can develop 
into a fully malignant tumour, is the best predictor for its evolving into a tumour. By studying the 
dynamics of a population described by a multi-type, population-size limited branching process in terms 
of the evolutionary formalism, we show some general principles regarding the probability of a resident 
population to being invaded by a mutant population in terms of the number of types present in the 
population and their resilience. We show that, although diversity in the mutant population poses a barrier 
for the emergence of the initial (benign) lession, under appropiate conditions, namely, the phenotypes 
in the mutant population being more resilient than those of the resident population, a more variable 
noeplastic population is more likely to be invaded by a more malignant one. Analysis of a model of 
gene regulatory networks suggest possible mechanisms giving rise to mutants with increased phenotypic 
diversity and robustness. We then go on to show how these results may help us to interpret some recent 
data regarding the evolution of Barrett's oesophagus into throat cancer. 

Author Summary 

Recent results by Maley and co-workers [T regarding progression of the Barrett's esophagus, a benign, 
neoplastic throat condition, into a malignant tumour have revealed that clonal diversity within diseased 
tissue is the best risk predictor outperforming usual genetic markers. So far, there is no detailed expla- 
nation as to why this is so. What is the role of diversity in the progression to cancer and what are the 
mechanisms involved? In this paper, we address these issues and put forward some generic mechanisms 
that may help to gain a better understanding of the empirical results obtained by Maley et al. Our ap- 
proach to this problem is in two steps. First, we propose a simple model of gene regulatory networks and 
analyse their behaviour upon gene silencing. We observe that, following gene silencing, both phenotypic 
diversity and robustness increase. We then proceed to analyse how these two effects affect the dynamics 
of the cell populations and find that more diverse pre-malignant lessions are likely to evolve into a more 
malignant form only if phenotypic robustness increases, hence suggesting that robustness is essential to 
the results put forward by Maley and coworkers. 

Introduction 

Complex diseases such as cancer pose an enormous scientific challenge. Cancer involves an extra level of 
complexity as its development involves evolutionary processes driven by Darwinian natural selection [2]. 
The traditional view of the emergence of cancer consists of a series of mutations with each of them 
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increasing the fitness of the mutated cells with respect to that of their normal counterparts, which are 
eventually taken over by succesive rounds of clonal expansion [3|[4j. Recent experimental analysis on the 
genetic landscape of tumours |5 , however, show an even more complex situation whereby tumours have 
been found to be much more genetically diverse than expected where low-frequency mutations are likely 
to play a leading part in the evolution of tumours. Furthermore, recent analysis of data from patients of 
Barrett's oesophagus, a neoplastic condition that may evolve into throat cancer, has found that the best 
predictor for its evolution into a fully malignant lession is diversity in the cell population, outscoring all 
of the usual genetic markers [T]. In particular, Maley et al. have shown that quantities normally used 
in ecology to quantify diversity, such as the Shannon entropy, Hs, and the number of clones detected 
in tissue samples extracted from patients of Barrett's oesophagus, are excellent predictors. According 
to [4], these findings bring into question the vision of tumour evolution as a series of succesive rounds 
of mutation/ clonal evolution to favour a model in which prcmalignant lessions sustain a number of co- 
existing cellular types. 

In order to explore some of the issues put forward by results of Maley et al. [l], in particular their 
results concerning the higher probability of pre-malignant lessions with larger diversity to evolve into 
tumours, we first analyse the robustness process of a simple model of gene regulatory networks. We 
focus on the study of the effects of gene silencing on both the number of stable phenotypes and their 
robustness [6}jll| . We will show using a simple model, similar to the ones used in previous related 
works [7]-|9{|12|, that such modification on the gene regulatory network leads to an increase on both the 
number of stable (robust) phenotypes and their ability to sustain further mutations (robustness). 

We then analyse how these changes induced by gene silencing have on the dynamics of the corre- 
sponding population, in particular on the ability of one mutant carrying such modifications to invade 
the resident cell population. To so doing, we consider a model of population dynamics consisting of a 



multi-type branching process with population-size limited proliferation probability 13 . Our aim is to 
put forward a number of general properties of multi-type populations which could shed some light on 
the mechanisms involved in the transition from neoplasm to malignant tumours. Our analysis is car- 



ried out in terms of the so-called evolutionary fornialim developed by Demetrius and coworkers 14 ■ 16 
This framework extends the thermodynamic formalism developed within ergodic theory and statistical 
physics [17| and it allows us to study the evolutionary behaviour of structured populations in terms of 
the evolutionary entropy^ which is a measure of the dynamical diversity of the population. 

The evolutionary formalism, whose main results concerning the present problem are summarised 
in Materials & Methods, has been recently applied to analyse the competition between resident and 
mutant populations modelled in terms of age-structured population models |16j and multi-type branching 
processes 13 . Demetrius et al. T5p6 have shown that the probability of a mutant structured population 
to invade the resident one is determined by the rate at which a population returns to its corresponding 
steady state after a random perturbation of the parameters that characterise the population dynamics. 



In turn, it has been proofed that such rate is given by the evolutionary entropy 15 . They have also 



shown that the evolutionary stability of a population, that is, its resilience to invasion by a mutant, is 



determined by the evolutionary entropy 16 



Our aim in this paper is to understand the dynamical mechanisms involved in the findings of Maley et 
al. [1]. In other words, we intend to find how the Shannon entropy, which is an index of diversity within 
the population, relates to the evolutionary entropy, which determines the evolutionary stability (or lack 
thereby) of the corresponding population. 
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Results 

Robustness of gene regulatory networks 

In order to model the dynamics of gene regulatory models, we consider a network with Nq nodes, each 
representing a gene. To each of the nodes we associate a state gi = ±l,i = 1,...,Nq, with state 
g = 1(— 1) corresponding to the gene being activated (inactivated). We generate a directed network, 
defined in terms of the corresponding adjacency matrix, according to a prescribed degree distribution 
(in the present case the in- and out-degrees are exponentially distributed). The links of this network 
correspond to the interactions between genes, which can be activatory or inhibitory. If a gene i activates 
gene j, the corresponding link, represented in a matrix W — (wij) whose entries are zero if there is no 
link between nodes i and j and non-zero otherwise, carries a weight Wij = 1. If, on the contrary, gene 
i inhibits gene j, the corresponding weight Wij = — 1. In our model, positive and negative weights are 
randomly assigned with probability p+ and 1 — p+, respectively. 

The dynamics of the gene regulatory netowrk is defined in terms of the above elements as follows. An 
initial condition gi(t = 0),i = 1, . . . , Nq and a matrix W are given. For each node and at each time step 
the following quantity is computed: Ii(t) = '^ji9ji^)- Then at each time step: 

r iif/,(i)>o 

g,{t + l) = \ -lif/,(t)<0 (1) 
[ g^{t) if /,(i) = 

which is ran until the system reaches a steady state. The phenotype, </>, is defined as the vector </> = 
[gi, . . . jgNc) steady state. This steady state may be a periodic solution. How we deal with this type 
of solution will be explained in more detail later. 

In our model, we consider two type of mutations. On the one hand, following [7|[9j[T2], we introduce 
random changes on how genes regulate each other (i.e. Wij — >■ —Wij with a certain probability). The 
second type of mutation we consider is a form of gene silencing where one particular gene is chosen (in 
principle, at random) and its state is fixed at gi = — 1 regardless of the state what remains of the network. 
This event could model both loss-of-function mutations j2j or an epigenetic silencing event [18 . 

In order to analyse the robustness properties of the corresponding phenotypes, 4>, which we define 
as the probability that a mutation of the first type according to the description in the paragraph above 
induces a phenotype switch (defined as a change in sign of at least one of the components of the vector 0) , 
we carry out a numerical experiment using the algorithm defined in the Materials & Methods section. Note 
that the population dynamics defined by this algorithm is completely neutral, since all the individuals 
yield offspring with the same probability and we have not imposed any viability conditions as assume e.g. 
in [t]. The corresponding results are shown in Fig. [l] We see that even in neutral conditions the diversity 
of phenotypes within the population decreases, while the extant types exhibit a greater resilience against 
mutations. In other words, the number of different phenotypes decreases with time but the surviving 



ones are such that their robustness increases, i.e. we observe the emergence of canalisation 18 19 within 
our simple model. 



These results can be understood in terms of the concept of the neutral network as introduced in 20|21 
The neutral network is a meta-network where each node corresponds to a different choice of W and two 
networks are connected if they differ only in one of the values of the corresponding entries Wij . It has 
been shown 20^ that phenotypic robustness can be understood in terms of the structure of such network: 
robust phenotypes are those whose subnetwork is highly interconnected, i.e. mutation events are likely 
to produce a network with the same phenotype. This means that robust phenotypes act like traps within 
this network making it harder for mutations to drive the system out of the corresponding phenotype. The 
same topological mechanism causes non-robust phenotypes to disappear at early stages in the evolution of 
the population. This is because mutations will drive them towards network configurations representing 
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the more stable phenotypes. Thus, we can quahtatively explain both why the number of phcnotypes 
decreases in time and also how the remaining phenotypes become more resilient. 

We now repeat the same simulatons with the difference that at a given point in time we randomly 
select one of the nodes, say node i, and introduce a gene silencing event, i.e. gi = —1 which we maintain 
fixed as the rest of the network evolves according to Eq. ([T]). The results are shown in Fig. [2j We 
observe that the induced mutation induces both an increase in the number of robust phenotypes and 
their robustness. These findings are consistent with experimental observations concerning liberation of 
cryptic genetic variation by inactivation of an evolutionary capacitor [6 as well as with recent modelling 
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Note that we do not explicitly deal with oscillatory asymptotic states of the gene regulatory network. 
In the present case, however, this is unimportant for the following reason. If from a fixed point a 
mutation generates an oscillation, this corresponds to a different phenotype and the phenotype switch 
will be correctly counted as such. If, previous to the mutation, the gene regulatory network is in an 
oscillatory state and the mutation does not produce a change in phenotype but simply, say, a phase shift, 
which does not have any bearing on the phenotype exhibited by the network, this will be counted as 
change in phenotype. This, however, makes our point on the increase in robustness by gene silencing 
even stronger, as what we are actually measuring is an lower bound of the phenotypic robustness. 

Now that we have established how gene silencing affects our model of gene regulatory networks, we 
move on to the analysis of effect of mutations on the dynamics at the level of cell populations. In 
particular, wc are interested in the competition between a resident (normal) population of cells and a 
population initiated by one cell that has undergone a gene silencing event. We investigate under which 
conditions the mutant will take over the resident. We start by introducing the main results of the 
mathematical formalism that we use to analyse this situation. 



Results for multi-type branching process with resource limitation. 

In order to advance further our discussion, we focus on a multi-type branching process with resource 
limtation. In the context of this paper, the different types in which the population divides correspond 
to the robust phenotypes. We further assume that all the (pheno)types yield offspring with the same 
probability, namely that of an individual in our original population model. For the sake of example, we 
will consider two different populations. The incumbent population is defined by an strategy whereby the 
population splits in d/ types. The average number of offspring to be produced by a given indvidual is 
given by: 



A = 2e-'^^ 



V 



-1 
V 



\ 



di-l 



(2) 



where Vij is the mutation probability per generation from type i to type j. For simplicity we will assume 
Vij — V Vi,j. Here N = Nj + Nm, i-e. the total population (both resident and mutant). In all the 
simulations shown in this paper the initial mutant population is Nm — 1- 

The mutant population adopts a strategy with a larger variety of available types {(Im > di different 
types) but which may have an increased robustness (smaller type mutation probability). Otherwise, the 
proliferation probability is the same as for the incumbent population. Thus: 



B = 2e-^^ 



dM-1 
dM-1 



(3) 
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InEq. ^, /9 is a measure of the increase in resilience of the corresponding phenotype. The parameter 
(j) corresponds to an increment in the growth rate of one of the phenotypes, which can be either positive 
or negative. The remaining phenotypes produce offspring at the same rate as the phenotypes of the 
incumbent population. 

With relation to the results of Maley et al. on the Barrett's esophagus |1 , the competition between the 
two populations in our model could represent the initial emergence of the neoplasm (mutant population) 
which overtakes the normal tissue (incumbent population). 

There are a number of results that we can obtain from numerical simulation of the corresponding 
branching process (see Materials & Methods for details). First, increased diversity {cIm > di) without 
increased resilience results in loss of the ability of the mutant to invade the resident population, this is 
because Pp decreases when cIm increases (see Fig. [s]). In other words, increasing the phenotypic diversity 
induces an error catastrophc-likc behaviour. Furthermore, increasing resilience (i.e. increasing p) rescues 
the more diverse mutant populations from the error catastrophe, so that they can now invade populations 
with less diversity, as shown in Fig. [3j Both these results can be explained in terms of the evolutionary 



formalism 13 . We can also observe (Fig. |4]) that Hs and Pp are negatively correlated: the larger 
the Shannon entropy of the mutant population the less likely fixation is. In other words, high-diversity 
mutant populations are less likely to invade and reach fixation. 

We consider now the stability and resilience to invasion of the mutants that have achieved fixation. 
Within the context of the problem we are dealing with, namely, the transition from neoplastic lesions to 
fully malignant tumours, this transition requires further mutations that perturbs the neoplastic cells and 
drives them into the malignant state. We model this type of mutations as perturbations in the parameters 
that determine the population dynamics. In particular, we consider that a second phenotype acquires a 
slight advantage in its proliferation probability quatified by a factor e'^. We assume that perturbations 
in the resilience of the phenotypes is represented by an additive term Ap. The mean- field dynamics (see 
Appendix) of this second mutant population is thus described by the matrix: 



B = 2e 



-fiN 



f e-il-. + p + Ap) e-^^ ... e^^-^ \ 

^^"^ eni-- + P + Ap) ... e-^^ 

L'—p—Ap u — p-\-ls.p 

(Im-I dM-l 

V ^ l-. + p-Ap) 



(4) 



We consider now essentially the same problem as above, namely, the likelyhood of a population whose 
(deterministic) dynamics is described by B to take over a resident population whose dynamics is given by 
B (Eq. ([3|). Using the methods described at length in our previous work 13 , the relative fitness s (see 



Materials & Methods for the precise definition of this quantity) of B with respect to B, which characterises 
the fixation probability of the former when competing with the latter, can be computed. The results are 
shown in Fig. [5] We observe that in the case in which the original mutation (that is, the one giving 
rise to the neoplastic lesion) increases the diversity of the population, i.e. djvf, without a concomitant 
increase in the phenotypic resilience, i.e. p = (see Fig. [5j red dashed line), the ability of the neoplasm 
for further evolution decreases with increasing cLm- In other words, the robustness of the neoplastic lesion 
increases with cIm , as the probability of it being invaded is a monotonically decreasing function of cIm ■ If, 
on the contrary, the original mutation, in addition to increasing the number of cellular types within the 
population, also increases their resilience, i.e. p > 0, the robustness of the resulting population decreases 
with increasing cIm, as the corresponding relative fitness, and therefore the fixation probability, is a 
monotinically increasing function of dm (see Fig. [s] black solid line). Numerical results corresponding 
to this situation are presented in Fig. |6] and show agreement with the theoretical predictions of the 
evolutionary formalism. 
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Discussion 

In this paper we have analysed the effects of gene silencing by means of gene silencing on a model of 
gene regulatory networks and we have studied how it affects phenotypic diversity and robustness. We 
have then investigated how these changes can modify the ability of a resident population to withstand 
invasion by mutants, and, last, we have discussed how our general results may shed some light on recent 
results regarding the evolution from neoplastic lesions to fully-malignant tumours. 

We have shown that, within our model, gene silencing leads to both an increase in phenotypic diversity 
and also enhances the robustness, i.e. in the populations ability to sustain gene mutations. Furthermore, 
we have done this in a completely neutral model with no viability conditions [jj. Our neutral model 
allows us to compute a lower bound on the robustness. 

We have then proceeded to analyse the consequences of the above observations for the competition 
between resident and mutant populations and on the ability of the mutant population to further evolve 
into more malignant variants. We have done this within the framework of multitype branching processe 
with resource limitation and by applying the evolutionary formalism. We have shown that increased di- 
versity, as measured by the corresponding Shannon entropy, correlates with decreased fixation likelyhood 
for the initial invasion (i.e. the emergence of the neoplasm, e.g. Barrett's oesophagus). If, in addtion 
to an increased number of phenotypes, there is a concomitant increase in the phenotypic resilience, i.e. 
p ^ 0, the likelyhood of invasion increases. However, even the latter case where invasion is more likely, 
the same general pattern of anti-correlation between Shannon entropy and fixation probability remains 
(see Figs. |3]and|4|. Furthermore, we have shown that both increased phenotypic resilience and increased 
diversity correlate with lower robustness of the population and therefore they are more likely to be further 
invaded. 

These general results on population dynamics can help us to make sense of some recent data regarding 
the progression of the Barrett's oesophagus into a fully malignant throat cancer [l]. First, less than 5% 
of the cases of Barrett's oesophagus actually evolved into cancer. This fact can be understood in terms 
of our results: lesions with increased diversity, which are the more likely ones to evolve into a tumour, 
are going to be less frequent, as their fixation probability is smaller. Furthermore, we predict that the 
fact that more diverse lesions are more likely to evolve into cancer implies that the event leading to the 
pre-malignant lession must also increase the phenotypic resilience, as that induces loss of robustness in 
more diverse populations. 

Materials and Methods 

Dynamics of the population of networks In order to analyse the robustness properties of the 
corresponding phenotypes, 0, which we define as the probability that a mutation of the first type according 
to the description in the paragraph above induces a phenotype switch (defined as a change in sign of 
at least one of the components of the vector (/>), we carry out a numerical experiment defined by the 
following algorithm: 

1. An initial population of A'o W-matrices, i.e. directed networks, are randomly generated as described 
above. Each of these matrices is assumed to define a particular individual within the population. 

2. The gene regulatory network dynamics described above is iterated for a number of time steps (fixed 
to a much larger value than typically needed to reach a steady state) for each of the Nq individuals. 
The vector (gi , . . . , gNa ) the end of iteration round is taken to be the phenotype corresponding 
to each individual within the population. 

3. The population is next binned according to the corresponding phenotype, with . being the 
number of individuals with a particular phenotype 4>j = {gi, i = 1, . . . , Nq} 
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4. Each individual produces two offspring with probabihty Po = a t^^W^ where /i is the carrying 
capacity 

5. Upon division, each regulatory interaction is mutated (wij — >■ —Wij) with probability p„ = rm/i?, 
where r„i is the mutation probability per generation per edge and E is the number of non-zero 
entries in the matrix Wij 

6. If no mutation occurs, two new individuals are incorporated in the next generation with the same 
Wij and 

7. If mutations happen, steps 1 to 3 are repeated and new individuals are created accordingly 



Summary of the evolutionary formalism. Here we present a brief summary of the main results and 



formulae of the evolutionary formalism. For a full account of the details we refer the reader to 14 ■ 16 



For specific results concerning the particular model we are analysing here, see 13 . 

Consider a population with dj different types whose dynamics is given by iteration according to 
u{t + 1) = A{t)u{t), where u{t) — {ni{t), . . . ,ndi{t))'^ is the vector consisting of the population sizes of 
each of the types and A = (aij) is the average number of offspring of type i produced per generation by 
individuals of type j. Note that, in general, the entries of this matrix depend on the total population 
Nj — J2iLi '^i- We will consider the long-time behaviour of the population, i.e. once it has settled in a 
steady state. We will further assume that A is an irreducible matrix. 

The long time behaviour of the population is determined by the dominant eigenvalue of the matrix 
A, which will denoted by Aq. Let u and v be the corresponding left and right eigenvectors, respectively: 
Au = AqU and vA = Aqv. We can now define an associated Markov matrix, P, such that: 

P^j = ^ (5) 

One of the main results consists in a variational principle, which states that the following relation 
between the growth rate r = log Ao, the evolutionary entropy, He, and the so-called proliferative potential, 



F, 14 



r = HE + F 

di 

He = ^ T^iPij log Pi j 

i,j = l 
d, 

F = ^ TTiPy logCij (6) 
i,i = l 

where tt^ is the stationary distribution corresponding to the Markov process defined by P. The requirement 
that A should be irreducible ensures the uniqueness of this distribution. 

A number of interesting results follow from Eq. (j6| and the corresponding variational principle, for 
example, a proof concerning large deviations states that the rate of relaxation to the steady state after 
a perturbation of the entries of the matrix A, specifically, a perturbation of the type A((5) = {a^^^), is 
positively correlated with the corresponding change in the evolutionary entropy He{5) — He{5 = 0), thus 



proving that this quantity is a measure of robustness of the population 15 . 

However, we are more interested in a (related) result concerning the evolutionary stability of an 
incumbent population against an invading mutant population when both compete for a common limited 



resource. By means of a diffusion approximation, Demetrius et al. 16 have showed that the fixation 
probability of the mutant, Pp is given by: 



where the total population is assumed to be constant, y is the initial proportion of mutants, and s is 
defined by: 



where (N) is the average stationary population. In the case of the branching process we consider in this 
paper, (N) is given by 2e~''^^^ = 1, i.e. IJ.{N) — log 2. The quantity Ar is defined as Ar = TM — rj, where 
the subindices / and M denote the corresponding quantities for the incumbent and mutant populations, 
respectively. The quantity Aa^ is defined as Aa^ = cr'^j — <t| where 13 : 



dHE{6) 



d6 



(9) 

(5=0 



The behaviour of Pp can be analysed in terms of the quantity s. If s > 0, Ppiv) is convex and 
Ppiu) > y- If, on the contrary, s < then Ppiv) is concave and Ppiv) < V- This means that if s < it 
is very likely that the mutant gets extinct whereas if s > it is very likely that the mutant is fixed, with 
the likelihood of invasion increasing as s increases. Moreover, the larger the value of s the more likely 
fixation is. Therefore s can be taken as a measure of fitness. 

Definition of the muitype branching processes To precisely define the branching processes in- 



volved in our models of population dynamics, we use the corresponding generating function formalism 22 
where a generating function for the probability distribution of the per individual offspring number is pre- 
scribed. In this article we consider three different populations (corresponding to the matrices A, B, and 
B). The set of generating functions corresponding to the population associated to the matrix A is: 

9,{z) = (1 - e-^^) + (1 - .)e-^^zf + J] ^e"-^z| (10) 

where i = 1, . . . , c?/ and z is a dj-dimensional vector with < z < 1. 

Similarly the generating function for the populations associated to B and B are, respectively: 

51 (z) = (1 - e-"^+'^) + p))e-''^+^z2 



E 



dM-l 



gm = (1 - e-^^) + (1 - (^^ - P))e-'^^zf + 

y ^^^e-A'^z2 fo, i > 2 (11) 
tl^dM-1 ^ " 



and 
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^^ z.-p-Ap ^_^^ for z < 2 
(l-e-^^) + (l-(j.-p-Ap))e-^^^f 



cJm - 1 



e-'^^zl for i > 2 



(12) 



For both Eqs. ( 11 1 and ( 12 ) where i = 1, . . . , (Im and z is now a c?M-diinensfonal vector with < z < 1. 



In our simulations, at each generation we go through all individuals and kill them with a probability 



given by the Oth-order term in the generating functions Eqs. (10 1, (11) and (12 1, thus passing zero 



individuals on to the next generation. Otherwise, two individuals are passed on to the next generation of 
the same type, say type i, with probability given by the coefficient of the yf- or z|-terms or of a different 
type, say type j, with probability equal to the coefficent of the yf- or z^-terms. 

Note that the entries matrices A, B, and B correspond to the branching ratios m^j = djgi{l) as defined 



in Eqs. (10), (11) and (12), respectively, and therefore they determine the mean-field dynamics. 
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Figure 1. Simulations of the evolution of robustness and phenotypic diversity in a model of gene 
regulatory networks (see text for details). The upper plot corresponds to the phenotypic sensitivity to 
mutations, i.e. the probability that a mutation induces a change in phcnotypc, which is an inverse 
measure of robustness: the more sensitive, the less robust the phenotype. The lower plot corresponds to 
the phenotypic abundance defined as the number of phenotypes that are actually present within the 
population at generation t divided by the total population. We observe that both quantities decrease 
with time. Parameter values: Nq = 10, /x = = 0.0002, = 1 
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Figure 2. Simulations of the evolution of robustness and phenotypic diversity in a model of gene 
regulatory networks (see text for details). The upper plot corresponds to the phenotypic sensitivity and 
the lower plot to the phenotypic abundance (see caption of Fig. [T]). In these simulations we have 
introduced a gene silencing event (see main text for details) at generation t = 250, which is released at 
t = 500. Parameter values: Ng = 10, u = K^^ = 0.0002, r™ = 1 




Figure 3. This figure sliows tlic fixation probability, Pp, as a function of (Im corresponding to = 0.01 
for jjL = 0.0001 with p = (squares) and p = 0.49 (circles) and n = 0.001 (squares) with 5 = 0. The inset 
shows the fixation probability for ip = 0.01 and fj, = 0.0001 as a function of the phenotypic resilience, p. 



14 



1.5 

o 

h 
c 

<u 

c 
o 

§ 1 

C/2 



0.5 



1 ^ r 



1 ' r 



h-«-H 



h-»- 



I • I 



\-m-\ 



I — ■ — I 



I • I 
I • I 



I — • — I 



I ■ 1 



I — • — I 



J I L 



_L 



J . ' I 



I r« H 



0.0025 0.005 0.0075 0.01 0.0125 0.015 



Figure 4. This plot shows the negative correlation between the average Shannon entropy of the 
mutant population and the corresponding fixation probability. Circles correspond to simulation results 
for p = and squares to p = 0.49. 




Figure 5. This plot shows the invasion fitness, s, as a function of cLm- The red dashed line and the 
black solid line correspond to p = and p = 0.49, respectively. Ap — 0.0075 in either case. 
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Figure 6. This plot shows the fixation probability of a population described by B over a B-population, 
s, as a function of cLm- Squares correspond to p = and circles to p = 0.49, respectively. Ap = 0.0075 
in either case. 



